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ABSTRACT 

'3CO(J=2-l) and C"*0(J=2-1) observations of the molecular cloud G285.90-l-4.53 (Cloud 16) in the Carina Flare supershell 
(GSH287+04-17) with the APEX telescope are presented. With an algorithm DENDROFIND we identify 51 fragments and com- 
pute their sizes and masses. We discuss their mass spectrum and interpret it as being the result of the shell fragmentation process 
described by the pressure assisted gravitational instability - PAGI. We conclude that the explanation of the clump mass function needs 
a combination of gravity with pressure external to the shell. 

Key words. Stars: formation; Stars: winds and outflows; ISM: molecules; ISM: structure; Galaxy: structure 



1. Introduction 

Studies of neutral a tomic hydrogen (HI) in the Milky Way 
dKalberla et alj2005h and in nearby galaxies (Walter et al. 2008) 
have discovered shells, fountains and chimneys as ubiquitous 
fine-scale features in the HI gas content of ga laxies (HeilesI 
119791; lEhlerova & PalousI l2005t iBagetakos et all lIoTTI) . In the 
majority of cases these structures may be explained by the feed- 
back from young and massive stars. Stellar radiation, winds and 
supernovae of OB associations and of super star clusters are 
the main driving force forming the ISM holes and expanding 
HI shells. In the HI searches, three kinds of cavity are dis- 
tinguished: complete holes surrounded by dense ISM shells, 
one-sided shells which have blown out of the galactic disk 
into the halo in only one direction, and chimneys, in which a 
shell has blown out on both sides, puncturing the galactic disk. 
These latter blown-out structures provide hot gas and metals 
to galactic haloes. The chimney walls formed from the swept- 
up gas may self-shield against UV photons, form molecules 
and cool to very low temperatures. With the help of the com- 
bined acti on of thermal instabilities ("Kovama & Inutsuk i2002t 
lAudit & H ennebelle 2005; Vazquez-Semadeni et al. 200 ^, hy- 
drodynamic instabilities ( Vishniac 1983; Rvu & Vishniac 1991") 
and gravitational insta bilities (Vishn iac 1983,; Elmegreen 1994; 
IWhitworfli" et al.llT994l) . the molecular layers can fragment into 
clumps, forming sites of secondary triggered star formati on. 

In a recent series o f papers (Dale et al. 2009; Wiinsc het alj 

I20T0I iDale et"ani20ril) . we investigated the gravitational frag- 
mentation of a shell expanding into a high-pressure but low- 
density medium confining the shell by thermal pressure but not 
by ram pressure. We observed that the ambient pressure was of 
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crucial importance in determining the size scales on which the 
shell fragmented. We now seek to confront our theoretical work 
with observations of fragments in the walls of the Carina Flare 
supershell. 

The Carina Flare supershell, GSH287H-04-17, on the near 
side of the Carina Arm is one of the closest examples of an ex- 
panding shell reaching high galactic latitudes in the Milky Way 
disk. It was discovered in ^^CO(J=1-0) with the NA NTEN tele- 
scope bv lFukui etall (ll999l) . lDawson et alj(l2008bl) used Parkes 
telescope HI observations to define the Carina Flare as an ellipti- 
cal shell approximately 230 x 360 pc centred ~ 250 pc above the 
Galactic Plane. They computed its HI mass as 7 ± 3 x 10^ M© and 
its molecular hydrogen mass as 2.0 +0.6 x 10^ M0, and estimated 
the shell's kinetic energy as 10^' erg. 

iDawson et aP (12008a) used the NANTEN millimetre tele- 
scope to study the molecular content in more detail. In 
'2CO(J=1-0) and '^CO(J=1-0) observations, they identified 
~220 clouds down to a completeness limit of ~ 60 M©. They de- 
rived fragment mass functions with slopes of 1.5 and found 

that all of their clumps had strongly subvirial masses. This latter 
result immediately raises the question of what instability caused 
the Carina Flare to fragment, since the action of the gravitational 
thin- shell instability would be expected to produce fragments 
with virial or super-virial masses. 

In a parsec-resolution study of the HI and CO distributions in 
GSH287-H04-17, Dawson et al. (2011b) distinguished between 
different origins of the molecular clouds in the Carina Flare 
walls. They identified the remnants of pre-existing clouds that 
have been engulfed by the hot medium in the expanding HI hole, 
as well as locations where molecular clouds may have condensed 
from the swept-up atomic medium collected into the shell walls. 
Of the latter, one of the most compelling examples is the molecu- 
lar cloud G285.90-h4.53 (hereafter Cloud 16, as in lDawson elaH 
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l2008ah . which lies approximately halfway up the western side of 
the shell, ~ 250 pc above the Galactic Plane in a region named 
"the approaching limb complex." 

The '2cO(J=l-0), i^CO(J=l-0) and C'^O(J=l-0) lines in 
this region have also been observed by the Mopra 22m telescope, 
providing detailed information on the distribution and properties 
of the molecular gas to a resolution of 0.6 pc, and confirming 
massive s tar formation activity in the southernmost regions of 
the cloud dPawson et al ]l2011ah . In this paper, we report the first 
results of our new high-resolution observations of '^CO(J=2-l) 
and C"*0(J=2-1) in the Carina Flare Cloud 16 with the APEX 
telescopqj With our signal-to-noise ratio, which is much better 
compared to the MOPRA data, and a spatial resolution of ~ 0.34 
pc, we identify 51 clumps. We discuss their sizes and masses, 
some of which are close to virial. We also discuss the possible 
origin of their mass distribution, which can not be created by a 
pure gravitational instability. 

2. Observations of the Carina Flare Cloud 16 

The observations were carried out with the Atacama Pathfinder 
Experiment (APEX) 12-m diameter antenna in December 2010 
and January 2011. Using the APEX-1 receiver of the Swedish 
Heterodyne Facility Instrument (SHeFI), the observations were 
done simultaneously at the frequencies of the '^CO(J=2-l) 
(vrest = 220.399 GHz) and Ci**0(J=2-l) (v,-es, = 219.560 GHz) 
lines. The Fast Fourier Transform Spectrometer backend was 
used with a bandwidth of 1 GHz divided into 2048 channels, cor- 
responding to a velocity resolution of 0.66 km s ' . At 220 GHz, 
the half-power beam width of the telescope is about 27" which at 
the adopted distance of the Carina Flare of 2.6 kpc corresponds 
to about 0.34 pc. The system temperatures were typically in the 
range 150-170 K; only a very small part of the data had signifi- 
cantly poorer quality with Tsys up to 650 K. 

To map the brightest parts of Cloud 16, we observed six 
10' X 10' rectangular areas - four located in an "L"-shaped lower 
region, and two in an upper region, separated by about 18' (see 
Figure [TJ. We used On-The-Fly (OTF) mode which allows the 
telescope to scan continuously along one direction while read- 
ing out data every 10" after 1 second integration intervals. The 
scanning pattern used was a "zigzag" through the regions which 
were tilted by 26° with respect to the RA axis. The second row 
was scanned in the opposite direction, and so on, with 10" spac- 
ing. This pixel size (10" x 10") corresponds to the spatial resolu- 
tion 0.18 pc X 0.18 pc at the distance of Carina Flare (2.6 kpc). 
Each map (except the top box of the upper region) was observed 
twice in two perpendicular directions. The neighbouring regions 
overlapped by one beam width. The duration of one map with 
OTF overheads was about 2.5 hours. The total observing time 
was then 22.2 hrs of which about 2 hrs were lost. The observing 
conditions were good to excellent with pwv less than 0.6 mm. 

The data were reduced according to the standard proce- 
dure using CLASS from the GILDAS software package devel- 
oped at IRAM. Bad scans were flagged and emission line-free 
channels were used to subtract in most cases first-order base- 
lines. Some spectra showed a wide residual double sinusoidal 
variation caused probably by a vibration of either the cover of 
the Cassegrain cabin entrance window or the cold h ead of the 
closed-cycle cooling machine dde Riicke et al.l 120061) . We used 

' This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institute fur Radioastronomie, the European Southern 
Observatory, and Onsala Space Observatory. 



the original velocity resolution of ~ 0.66 kms The corrected 
antenna temperatures, T\, provided by the APEX calibration 
pipeline (Dumke & Mac-Auliffe 2010), were converted to main- 
beam brightness temperature by T\, - T^/ri^b, using the main 
beam efficiency 7],^^ = 0.75 at 220 GHz. Finally, FITS cubes 
were produced by CLASS for individual observed regions and 
stitched together manually using a Python routine. RMS noise 
levels of better than 0.2 K per 0.66 kms ' channels were ob- 
tained. 

The high angular resolution of APEX allows us in principle 
to observe objects of mass ~ 0.6 Mq, well into the regime of stel- 
lar masses. In Figure[T] we compare APEX and NANTEN veloc- 
ity channel maps to illustrate the improvement in resolution af- 
forded by APEX. Our observations show clearly that Cloud 16 is 
composed of irregularly shaped clumps connected by filaments, 
and that structure exists on scales considerably smaller than were 
accessible with NANTEN. 

3. Identification of Clumps 

The structure of the mo lecular gas in galaxies is hierarchical. 
iRosolowskv et al.l ( l2008l) showed that it is useful to display its 
essential features as dendrograms - graphical representations of 
the topology of the isosurfaces as a function of contour level. 
We have developed a simple algorithm called DENDROFIND 
which creates dendrograms as a natural by-product of the clump- 
finding process . The algorithm is similar to the well-known 
CLUMPFIND (IWilhams et al.llI994h but its results are less de- 
pendent on technical parameters (for instance, the brightness 
temperature difference between contours). Another advantage 
comparing to CLUMPFIND is that DENDROFIND always pro- 
duces contiguous clumps (CLUMPFIND may lead to discon- 
tiguous structures identified as a single clump in extreme cases). 

The DENDROFIND algorithm operates on a position- 
position-velocity datacube of the brightness temperature T\,. The 
main loop gradually decreases the temperature level Tl start- 
ing from the maximum brightness temperature in the datacube, 
Tmax, going down to Tcutoff set by the user. In each itera- 
tion, pixels with Tb > Tl are assigned to existing clumps, or 
new clumps are created and the pixels are assigned to them. 
Furthermore, information about connections among clumps on 
level Tl is evaluated and recorded. This information is neces- 
sary (and sufficient) to plot the dendrogram at the end of the 
algorithm. Each clump created by the algorithm has to fulfil the 
following three conditions: (i) it is a contiguous area in the PPV 
space; (ii) it consists of at least Npxmin pixels (Npxmin is an- 
other user defined parameter); and (iii) the difference between 
the clump peak temperature and the temperature at which the 
clump connects to other clumps has to be at least dTleaf . This 
parameter effectively controls whether a clump should be split 
into several objects, or just regarded as a single object with un- 
resolved substructure. 

In summary, the algorithm is controlled by four parameters: 
Nlevels (number of temperature levels), Npxmin (minimum 
number of pixels of a clump), dTleaf (minimum difference be- 
tween the clump peak temperature and the temperature at which 
the clump connects to other clumps) and Tcutoff (minimum 
Tb considered for assignment to clumps). It is always possible to 
find Nlevels such that using higher values would lead to identi- 
cal results. The DENDROFIND algorithm is described in detail 
in Appendix lAl 

In this work, DENDROFIND has been applied to the 
'^CO(J=2-l) and C'**0(J=2-1) data of the observed part of the 
Carina Flare Cloud 16. The parameters were Nlevels = 1000, 
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APEX: C^''0(J=2-1) APEX: "C0(J=2-1) APEX: ^^C0(J=2-1 ) NANTEN: "CO(J=1-0) NANTEN: ^^CO{J=1-0) rk- km 

(convolved to NANTEN r9SD[ution) [r\ Mri b J 
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I I I I I 

Fig. 1. Velocity integrated C'**0(J=2-1) and '^CO(J=2-l) emission observed with APEX (left two panels) versus '^CO(J=1-0) and 
^^CO(J=1-0) emission observed with NANTEN (right two panels). The middle panel shows the APEX '^CO(J=2-l) convolved to 
the resolution of NANTEN. The circle in the upper right corner corresponds to the beam FWHP (full width at half power). 




16D1-35 16C1 16B1-3 



Fig. 2. Dendrogram of the clumps in the '^CO(J=2-l) emission 
in the lower region. Colors of lines representing clumps are the 
same as colors of regions showing clumps in Figuresl4land lB.H - 
IB.3I Solid lines and filled regions (in Figures |4] and IB.H4B.3l l 
show clumps 16D1-18, 16C1 and 16B1-3, dashed lines and dot- 
ted regions (in Figuresl4land lB.ll4B.3l ) show clumps 16D 19-35 
(see beginning of ^for the nomenclature definition). 




16E1-9 16F1-3 



Fig. 3. Dendrogram of the clumps in the '''CO(J=2-l) emission 
in the upper region. Colors of lines representing clumps are the 
same as colors of regions showing clumps in Figures lBT4l4B.5l 



Npxmin = 5 and dlleaf = Tcutoff — 3crj^oise where iXnoise is 
the standard deviation of the noise in the data cube. The result- 
ing dendrograms are seen in Figures Hand [3 In the '^CO data, 
we identified 5 1 clumps, which are listed in Tables[T]and[3] Some 
of these clumps were also found in the C'^O data (see Tables |2] 
and nil. Figure m shows an example radial velocity channel map 
with the identified clumps. All radial velocity channel maps with 
clumps are shown in the online appendix iBl 



4. Clump Sizes, Velocity Dispersions and Masses 

Tables [T]-2] list properties of clumps identified in '^CO(J=2-l) 
and C'^0(J=2-1) APEX observatio ns. "Label" is the c lump- 
identifier based on the list given by Dawson et al. (2008a). The 
first nu mber and letter is identical to the list of Dawson et al] 
(l2008ab . e.g. " 16D", then a number denoting the substructure 
which we resolve is appended. "FWHM j" is the full width at half 
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Fig. 4. Clumps identified in the lower region (left two panels) and corresponding brightness temperature (right two panels) in the 
velocity channel v - 23.4 km s ' . Filled regions show clumps 1-18, dotted regions show clumps 19-35. Colors of clumps correspond 
to colors of dendrogram leaves in Figure |2] 



maximum of the clump in the direction of coordinate x, where 
X e {I, b, V = vlsr)- FWHM V is given by the formula 



FWHM^ = 2 V21n2 



Ax 



2V21n2 



(1) 



where the second term under the square root represents the pixel 
resolution term with Ax being the pixel width in the x coordinate. 
This term ensures the FWHM is not smaller than the pixel size. 
If the pixel size goes to zero the FWHM approaches a Gaussian 
FWHM. cr^ is the weighted variance, with the brightness tem- 
perature Tb as the weight: 



0-, 



_ — _ -2 _ X Tb _ ^—1 XI b 



xTb 



where T,r 



(2) 



(3) 



px 



The summation in the equations above is over all a given clump's 
pixels. R is the clump equivalent radius defined as 



(4) 



where is the clump projected area in the sky (/^7-plane). Tpeak 
is the brightness temperature of the brightest pixel of the clump, 
^peak, ^peak are its J2000 galactic coordinates and Vpeak is its LSR 
radial velocity. The next column, m, is the total clump mass com- 
puted as 



X 



where X = 0.735 



(5) 



and miij is the total clump mass in the form of molecular hydro- 
gen computed from the formula 



(6) 



px 



where fiun is the mass of one hydrogen molecule, Nh^ is the 
column density of molecular hydrogen in a given pixel and 
Apx is the pixel area. The column density A^H2 is derived from 
the column densities of the '^C O and C'^O molecules as - 
suming the standard abundances dBinnev & Merrifieldl [T998h : 
A^H, = 6.5 X lO-'^A^nco, A^Hj = 5 X 10^A?c>*o- These values are 
in agreement with recent measurements and models; see, e.g. , 
FSonnentruc keretaLl (l2007l): IVisser et alj (l2009h : IWilso nI (1999|) 
and the citations therein. Abundances, however, vary for differ- 
ent regions/clouds and depend also on the physical conditions 
inside individual clumps. Although this variation is a very funda- 
mental problem, a comprehensive discussion is beyond the scope 
of this paper 

The column densities of both '^CO and C'**0 are, in the 
case of 7 = 2-1 lin es, computed from the same formula 
(iRohlfs & Wilsonl[T996l 14.54) 



1-2 



1.51 X 10' 



1 _ e-io.6/7-„ j 



T dv 



(7) 



This formula assumes local thermodynamic equilibrium charac- 
terised by the excitation temperature Tex, which we assume to be 
25 K. This value is slightly higher than that usually assumed in 
such analyses (between 10 K and 20 K) and the value we use in 
]to analyze the fragmentation process (15 K). The reason for 
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this assumption is that the brightest pixel in our data is over 18 K 
(see Tpeak of clump 16D1), which puts the lower limit to the ex- 
citation temperature somewhere above 20 K. This is so because 
the excitation temperature has to be higher than the brightness 
temperature, assuming local thermodynamic equilibrium. Note 
that the column density, and thus the mass, is only higher by 
about 5 % for most clumps for an excitation temperature of 25 K 
compared to 20 K. In our case of discrete pixels, dv corresponds 
to Av, where Av is the pixel width in coordinate v (radial veloc- 
ity channel width, Av ^ 0.66 kms '). t is the optical depth of 
'3CO(J=2-l) or C''*0(J=2-1) given by the formula 



Label ^ 


r 

peak 


»,<") 


„,|ia) 


,„(13)_,„(18) 




K 


Mo 


Mo 


'"mean 


Lower region 


16D1 18.3 


3.2 


367.2 


296.7 


+0.21 


16D2 14.7 


2.5 


194.5 


168.0 


+0.14 


16D7 10.0 


2.5 


125.6 


139.1 


-0.10 


16C1 8.2 


1.4 


23.4 


22.6 


+0.03 


Upper region 










16E1 4.0 


1.5 


16.5 


36.9 


-0.76 


16E2 4.3 


1.6 


10.0 


19.8 


-0.66 



T = - In 



1 - 



T 

Io!6 



1 



exp(10.6/rex) - 1 



-0.02 



(8) 



Table 5. Comparison of C'^O clump masses derived from C'^O 
('18') and '''CO ('13') brightness temperature maps. 



which we derived from formula 14.46 of iRohlfs & WilsonI 
(Il996h . 

In addition to uncertainties caused by departures from LTE, 
the accuracy of the clump mass m depends on the data noise 
o-j. In our datasets, ctt is between 0.18 K and 0.20 K, which 
corresponds to an error in the determination of the pixel mass 
Ampx = 0.04 Mo for the '^^CO(J=2-l) line and Awpx = 0.3 Mq g DisCUSSion 



order Taylor expansion assuming the clump's area and FWHM,, 
errors of one pixel. The last column in Tables[T}Hlgives the ratio 
of the total clump mass derived from the brightness temperature 
and the clump virial mass. 



for the C"^0(J=2-1) line. Assuming the noise is Gaussian and 
neglecting the correlation of pixels within the beam, the error 
in determination of the clump mass (given in 11th column of 
Tables [Ba is 



Am = yA^pxAmpx 



(9) 



where A^px is the number of pixels in the clump. Since each 
clump has to have at least 5 pixels, the minimum clump mass 
is 0.6 Mq and 4.5 Mq for '^CO and C'^O data, respectively. 

To estimate the overall error in the mass determination, we 
compare our APEX '^^CO(2-l) observation with the NANTEN 
'^CO(l-O) observation by convolving the APEX data to the res- 
olution of NANTEN (see the middle panel of Figure [TJ. Both 
distributions of the brightness temperature show a very good 
agreement. The total mass detected in both lower and upper re- 
gions in the convolved APEX data is ~ 2.8 x 10^ Mq, which 
differs by only 15 % from the total mass detected in the same 
two regions in the NANTEN data (~ 2.4 x 10^ Mq). These val- 
ues also correspond well to the total mass in the identified '^^CO 
clumps (see Tables [T] and [3]) which is ~ 2.4 x 10^ Mq. 

As another test of reliability of the clump mass derivation, 
we compare masses obtained from the two lines observed by 
APEX. We select areas identified as clumps in the C"*0 data 
(Tables |2] and |4| and calculate masses contained in them from 
both C'^O and '^CO datasets. The results of the comparison are 
shown in Table |5] There are four C'^O clumps in the lower re- 
gion and two in the upper region. In the lower region, the mean 
difference between C'^O and '^CO masses is about 12 %, while 
in the upper region, the differences for the two clumps are about 
70 %. Part of these differences can be due to varying isotope 
abundances, different excitation temperatures in the lower and 
upper regions, and departures from LTE. From this comparison 
and from the comparison with the NANTEN data, we conclude 
that systematic errors in our determination of clump masses are 
probably less than a factor of two. 

The 12th column in Tables [T] ffl Wyir, is the estimated clump 
virial mass, given by the formula ( MacLaren et alJ[T988h 



As seen in Tables ri]^ the observed mas ses of individual clumps 
are, compared to iDawson et al.l (l2008 ah. significantly closer to 
virial, but they are still lower than the virial estimates. This dif- 
ference could mean that besides the Jeans gravitational instabil- 
ity, there are other forces helping gravity during the fragmenta- 
tion process. Another explanation may be that we are still limited 
in angular and velocity resolution, so that some of our clumps 
may be resolvable still further into smaller denser supervirial ob- 
jects. In particular, the relatively coarse velocity resolution is the 
main source of uncertainty in determining virial masses. This 
may be checked with future even higher resolution observations 
with ALMA interferometer. 



5.1. Clump mass function - CMF 

We have determined the mass function of clumps (CMF) iden- 
tified in '^CO(J=2-l) observations of both regions. Clumps ly- 
ing at the border of the observed area having Tpeak - T^boidei > 
dTleaf, where Tborder is the maximum brightness temperature of 
the clump at the border of the observed area, are excluded from 
the CMF, because their parts lying outside of the observed region 
may substantially contribute to their masses. There are 7 such 
clumps and they are marked by asterisks in Tables[T}Hl The CMF 
of the remaining 44 clumps is shown by Figure |5] The maxi- 
mum of the distribution is around mcMF ~ 10 Mq. The slope of 
the CMF is compared to the power law with the "Salpeter" slope 
-2.35 (dotted line) and a power law with slope - 1 .7 (dash-dotted 
line) found bv lKramer et al.l (11998 ') for clumps in several molec- 
ular clouds. The dashed Une shows the mass spectrum obtained 
from the Pressure Assisted Gravitational Instability (PAGI; see 



mvir = 190 ■ ■ FWHM,2 



(10) 



The formula above assumes spherical clumps with a density pro- 
file oc 1/r. The virial mass error Am,,,, is computed from the first 



The slope of the CMF agrees well with the - 1 .7 power law. 
The agreement with the PAGI mass spectrum slope (around -2) 
is a bit worse but still acceptably good. One reason why the 
observed mass spectrum is slightly less steep than the predic- 
tion of PAGI may be that the three most massive clumps include 
very dense regions that are not resolved completely. With higher 
resolution, these three clumps could perhaps split into smaller 
objects. At the low-mass end, the observed CMF includes four 
clumps with masses substantially lower than mcMF, whereas the 
PAGI mass function has a low-mass cutoff". One reason for this 
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Label 


FWHM_, 


FWHM,. 


FWHM,. 


R_ 










m 


Am 


'"vir 




m 


PC 


PC 


km s"' 


PC 


K 


de.E 


licR 


km S-' 


Mo 


Mo 


mT 


"ST 


mvir 


16D1 


0.59 


0.37 


1.98 


0.4 


18.3 


285:54:51.7 


4:31:16.6 


-22.75 


144.6 


0.33 


330.7 


221.4 


0.437 


16D2 


0.51 


0.51 


1.91 


0.5 


17.4 


285:54:37.4 


4:31:02.3 


-22.75 


119.6 


0.37 


315.2 


218.9 


0.379 


16D3 


1.08 


0.90 


2.07 


1.1 


15.5 


285:55:05.9 


4:31:45.2 


-23.41 


312.2 


0.79 


895.8 


574.0 


0.348 



Table 1. CO(J=2-l) clump properties - lower region. The meaning of the columns and the means by which the values were 
calculated are explained in ^ Clumps marked by the asterisk at their label are excluded from the CMF, because they lie at the 
border of the observed area. Full version of this table is available in the online appendix - see Table IB. ll 
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A/)( 
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PC 
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Mo 




"mT 


"Ivi,- 


16D1 


1.04 


0.79 


1.73 


0.8 


3.2 


285:54:51.7 


4:31:16.6 


-22.75 


296.7 


3.55 


481.9 


368.9 


0.616 


16D2 


0.59 


0.80 


1.51 


0.6 


2.5 


285:53:54.6 


4:30:19.3 


-23.41 


168.0 


2.77 


277.8 


244.4 


0.605 


16D7 


0.64 


1.29 


1.01 


0.7 


2.5 


285:54:08.4 


4:34:07.7 


-22.75 


139.1 


2.72 


138.8 


181.8 


1.002 


16C1 


0.41 


0.58 


0.66 


0.4 


1.4 


285:51:47.0 


4:23:52.9 


-25.40 


22.6 


1.12 


31.8 


63.6 


0.711 



Table 2. C'^0(J=2-1) clump properties - lower region. The meaning of the columns and the means by which the values were 
calculated are explained in ^ 
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16E1 


0.94 


1.66 


1.20 


1.2 


4.3 


286:00:20.5 


5:14:31.2 


-24.74 


121.0 


0.69 


335.8 


372.1 


0.360 


16E2 


1.32 


0.99 


0.88 


1.0 


4.2 


285:59:51.5 


5:13:20.1 


-24.08 


57.5 


0.47 


153.3 


231.1 


0.375 


16E3 


1.32 


1.29 


1.13 


1.1 


3.7 


285:59:10.6 


5:19:45.3 


-24.08 


56.4 


0.54 


270.0 


317.8 


0.209 



Table 3. '^CO(J=2-l) clump properties - upper region. The meaning of the columns and the means by which the values were 
calculated are explained in ^ Clumps marked by the asterisk at their label are excluded from the CMF, because they lie at the 
border of the observed area. Full version of this table is available in the online appendix - see Table IB. 2l 



may be incompleteness. We are unable to calculate the complete- 
ness limit without making further assumptions about the internal 
density and velocity structure of our clumps. Since our mass res- 
olution of ~ 0.6 Mq is, however, substantially lower than the 
maximum of the CMF at ~ 10 Mq, it is more probable that there 
really exists a small number of low-mass clumps that cannot be 
explained by the PAGI dispersion relation. They may be result 
of non-linear processes or non-gravitational fragmentation. 



5.2. Spatial distribution of clumps 

We estimate the typical separation of clumps, c/mst, from the av- 
erage length of the edges of the minimum spanning tree (MST, 
see Figure|6]l. Each clump is represented by the point at the po- 
sition of the brightest pixel of the clump (/peak, ^peak)- In each of 
the two observed regions, all the points are connected with lines 
- tree edges - in such a way that the sum of all edges lengths is 
minimal. Taking into account the distance to the Carina Flare, 
2.6 kpc, the average edge length in the lb-plane corresponds 
to ~ 1.4 pc with a standard deviation of ~ 1 pc. Since we do 
not have information about the positions of clumps in the ra- 
dial direction, we assume that their separations in this direction 
are similar to those in the other two directions, and multiply the 
average MST edge length in the lb-plane by factor V3/2. We 
obtain 

fl'MST = (1.7±1.2)pc. (11) 



5.3. Clump relative velocities 

Figure |7] shows the radial velocity distribution of clumps found 
in the ^^CO data. All clumps are distributed across only seven 
velocity channels, i.e. the maximum relative radial velocity of 



\-1.7 



\ \ 
\\ 



-2.35 



\ \ 



1 1.5 2 



Fig. 5. Mass spectrum of fragments identified in our '^CO(J=2- 
1) observations (solid line). Following the prescription of 
iMafz Apellaniz & Ubedal (l2005l) to minimize errors in his- 
tograms with a small number of objects, we set bin widths so 
that each bin includes approximately the same number of objects 
(six or seven). The dashed line shows the mass spectrum cal- 
culated from the PAGI dispersion relation (see equation \T7\ for 
I = gxlO-'^gcm-^, T = 15KandP,^^ = 1.4x lO^'^ dynecm^^ 
(parameters determined in ^5.4l i. Dash-dotted and dotted lines 
show power law with slopes -1.7 and -2.35, respectively. 



clumps is 4.6 kms"'. In some regions, the velocity varies by 
~ 2 kms ' on scales smaller than several pc of the projected 
distance, and there is no obvious global gradient in any direc- 
tion. This velocity variation is much higher than would be ex- 
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0.5 


1.5 


286:00:20.6 
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-24.74 


36.9 


1.53 


63.9 


102.1 


0.577 


16E2 


0.46 


0.37 


0.82 


0.3 


1.6 


286:00:20.5 


5:14:17.0 


-24.74 


19.8 


1.04 


41.4 


66.7 


0.479 



Table 4. C'**0(J=2-1) clump properties - upper region. The meaning of the columns and the means by which the values were 
calculated are explained in ^ 
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286.1 286 285.9 285.8 

I 



Fig. 6. Left: minimum spanning trees showing projected dis- 
tances between clumps. The two trees are constructed separately 
in the upper and in the lower observed region. Each tree node 
(denoted by the x symbol) is at the position of the pixel with the 
highest brightness temperatures in a given clump (Tpeak)- Right: 
outlines of the '^CO clumps projected onto the lb-plane. Colors 
and line types correspond to those in FigureslZjffland lBTljB.SI 



pected from the stretching of the shell surface due to the shell 
expansion. For instance, on a scale of 5 pc the relative velocity 
of clumps due to shell stretching is ~ 0.3 kms"' assuming the 
shell radius and expansion velocity to be 150 pc and 10 kms"', 
respectively, (Fukui et al. 1999). It may be, however, that the 
clump relative velocities are caused by their mutual gravitational 
attraction. For instance, if we estimate the dynamical velocity for 
Cloud 16D, we get Vdy„,i6D = (GMi6d/^16d)'^^ ~ 1.2 kms"' us- 
ing Mi6D ~ 1.8 X lO"* Mo and /?i6d ~ 5 pc. Another source of 
clump-clump relative velocities may be local processes like star 
formation and related feedback. Our spatial resolution, however, 
does not allow us to see these small scale processes. 



5.4. Origin of clumps: RAG/ 

The Pressure Assisted Gravita ti onal Instabilit y (PAG I) is 



described in Dale et al] (l2009h . IWiinsch et al] (I20I0}) . and 



iDale et al.l (l20I Ih . For PAGI. it is assumed that a thick shell ex- 
pands into a relatively rarefied ambient medium with non-zero 
pressure, P^^.^, and that the same pressure acts also on the in- 
ner wall of the shell. The shell fragmentation is controlled by 
the interplay between the gravitational force in the plane of the 
shell, the thermal pressure inside the shell and the external pres- 




V [km s"^] ^■'^ 



286.1 286 285.9 285.8 




-30 -28 -26 -24 -22 -20 -18 
V [km s"'[ 




Fig. 7. Distribution of the velocity in the '^CO APEX data. Top 
left: average velocity taken over all velocity channels weighted 
by Tb. Top right: integrated along the I coordinate shown in 
the ftv-plane. Bottom right: integrated along the b coordinate 
shown in the Zv-plane. Bottom left: the line profile integrated 
over both spatial coordinates. 

sure, Pext- PAGI can be described by the following dispersion 
relation: 



2,2 



lOP c^l 



5cH- 



2,2 



(12) 



where oj is the growth rate of a perturbation with the dimen- 
sionless wavenumber I. R, 2 and are the shell radius, surface 
density and internal speed of sound. The geometrical factor F re- 
sults from the approximation of a fragment by a uniform oblate 
spheroid and is given by 



where 



cos-'yS /3 

(1 -;62)3/2 " 1 _yj2 



r 



2c2 



I 



2P„,, + nGI? nR 



(13) 



(14) 



is the fragment aspect ratio (see lWiinsch et al.ll2010l for details). 
In Equation (fT2l) . we do not show other terms dependent on the 
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expansion velocity of the shell V, since the contribution of these 
terms is small compared to those from terms describing effects of 
the self-gravity, the external pressure, and the internal pressure, 
i. e., the first, second and third term of (fTSI) respectively. We also 
assume that the average molecular weight ju is 2.35 a.m.u. 

This dispersion relation produces a mass spectrum of frag- 
ments as 



diV 
dm 



oc o) X m 



(15) 



where m is the mass of an individual fragment and dA^ is the 
number of fragments in the mass interval (m, m + d m). A more 
detailed derivation of the mass spectrum is given in iDale et all 
(l20Tlh . 

The PAGI mass spectrum is in general a function of the frag- 
ment mass, m, the surface density of the shell, S, the external 
pressure, /"ext' ^^^11 temperature, T. Solid isolines in 

Figure |8] show the maximum of the mass spectrum m p^ai; as a 
funct i on of £ and f cvt- In h ne with previous works dFukui et alj 
ll999HDawson et al.i2008ah . we assume that the average temper- 
ature in clumps (and hence in the shell) is T = 15 K. This tem- 
perature is lower than the minimum we obtain from '^CO(J=2- 
1) observations in the densest regions (above 20 K) but it is prob- 
able that the temperature in these regions is higher than aver- 
age. In the derivation of PAGI, it is assumed that the mass of 
each fragment originates from a circular area with a diameter 
d = InR/l. The typical separation among neighbouring frag- 
ments is (ipeak - (4OTpeak/^S)'^^, shown by dashed isolines in 
Figure [8] An intersection between an mpeak isoline and a lipeak 
isoline gives 2 and P^^t of shell which forms fragments with 
typical mass mpeak and typical separation t/peak- 

Note that the intersection exists only for certain combina- 
tions of mpeak and lipeak- For instance, the mpeak - 3 M© isoline 
never intersects with the lipeak = 1 pc Such fragments are in- 
consistent with the PAGI theory because their surface densities 
are too low for fragmentation. Similarly, the mpeak - 30 Mq iso- 
line also never intersects with the t/peak = 1 pc but in this case, 
the surface density of such fragments is too high and they would 
break up into smaller objects. Therefore, this analysis provides 
quite a strong test if fragments with given masses and separa- 
tions can be created by PAGI or not. 

Taking the properties of clumps obtained from the '^^CO 
APEX observations (mpeak = mcMF - 10 and Jpeak - 
dusT =1.7 pc), we see that the corresponding isolines intersect 
at 2 = 9 X 10""* gcm"^ and P^^^ = 1.4 x 10"'^ dynecm"^ (the 
isolines are shown by the thick lines and resulting E and P^^^ by 
thin dotted lines in Figure [8]l. Considering the uncertainties in 
determinations of mcMF and c/mst, the clump properties would 
be also consistent with Pj,^^ > 1.4 x 10 dynecm"^ where the 
isolines are parallel, close to each other 

The obtained value of surface density is in an excellent 
agreement with the average surface density of the Carina Flare 
supershell determined in previous works (Z - MjiAnF?) ~ 
7 X lO"'* gcm"^ with the total mass of the shell M ~ lO*" Mq 
and its radius P ~ 150 pc; see iPawson et al]l2008bl) . The de- 
termined external pressure is somewhat higher that the thermal 
pressure thought to be present at this latitude (corresponding to 
z ~ 250 pc) above the galactic plane. The HI disc, however, 
cannot be supported against gravity by thermal pressure alone 
and other means of supp ort have been suggested, e.g., turbulence 
dPickev & Lockman|[l990i) . 



6. Conclusions 

We find that fragmentation in Cloud 16 of the Carina Flare 
has proceeded to much smaller scales than were visible in the 
NANTEN observations of Dawson et al. (2008a). It is probable, 
however, that there are structures present in this region that our 
APEX observations cannot resolve. The superior angular and ve- 
locity resolution of ALMA may allow the limits of the fragmen- 
tation process to be observed and accretion flows feeding clumps 
to be disentangled. 

Our th eoretical work dPale et al.1120091: IWiinsch et alj|2010l: 
iDale et al.>201 1 , PAGI) predicts that, in the early linear stages of 
shell fragmentation, the mass and the wavelength at which the 
shell preferentially fragments depends strongly on the confining 
pressure. Here we show that the CMF and the typical clump sep- 
aration in the observed part of the Cloud 16 are consistent with 
the fragmentation caused by PAGI. The surface density of the 
fragmenting layer, determined from the observed clump proper- 
ties, is E = 9 X 10""* gcm"^, a value in good agreement with the 
estimated average surface density of the Carina Flare supershell. 
The clump properties also suggest that the pressure confining 
the shell should be ~ 1.4 x 10 dynecm"^ or higher, we are 
unable, however, to exclude the possibility that the shell is in 
more advanced stages of the fragmentation for which the PAGI 
dispersion relation is not accurate. 

Since the extent of the Carina Flare above the Galactic Plane 
is several scale heights, the pressure encountered by the top part 
of the shell may be several times smaller than that confining the 
sides. We intend to perform observations of the molecular cloud 
G287.60h-8.00 (Cloud 74 of lDawson et"air2008 a). located at the 
greatest z-extent of the supershell, to see if the effect of the lower 
pressure environment on fragment size can be detected there. 

Acknowledgements. We thank the anonymous referee and Steve Shore for 
many suggestions which significantly improved the paper. We thank Michael 
Dumke for his assistance during APEX observations. This study was sup- 
ported by the Institutional Research Plan AVOZ 10030501 of the Academy of 
Sciences of the Czech Republic and project LC06014 Centre for Theoretical 
Astrophysics of the Ministiy of Education, Youth and Sports of the Czech 
Republic. VS acknowledges support from Doctoral grant of the Czech Science 
Foundation No. 205/09/H033 and a grant of the Charles University in Prague 
No. SVV 261301. 



References 

Audit, E. & Hennebelle, P. 2005, A&A, 433, 1 

Bagetakos, I., Brinks, E., Walter, E, et al. 2011, AJ, 141, 23 

Binney, J. & Merrifield, M. 1998, Galactic Astronomy 

Dale, J. E., Wiinsch, R., Smith, R. J., Whitworth, A., & Palous, J. 201 1, MNRAS, 
411,2230 

Dale, J. E., Wiinsch, R., Whitworth, A., & Palous, J. 2009, MNRAS, 398, 1537 
Dawson, J. R., Kawamura, A., Mizuno, N., Onishi, T, & Fukui, Y. 2008a, PASJ, 
60, 1297 

Dawson, J. R., McClure-Grifliths, N. M., Dickey, J. M., & Fukui, Y. 201 la, ApJ, 
741, 85 

Dawson, J. R., McClure-Griffiths, N. M., Kawamura, A., et al. 2011b, ApJ, 728, 
127 

Dawson, J. R., Mizuno, N., Onishi, T, McClure-Grifiiths, N. M., & Fukui, Y. 

2008b, MNRAS, 387, 31 
de Rijcke, S., Buyle, R, Cannon, J., et al. 2006, A&A, 454, LI 1 1 
Dickey, J. M. & Lockman, E J. 1990, ARA&A, 28, 215 

Dumke, M. & Mac-Aulift'e, F. 2010, in Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, Vol. 7737, Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series 

Ehlerova, S. & Palous, J. 2005, A&A, 437, 101 

Elmegreen, B. G. 1994, ApJ, 427, 384 

Fukui, Y, Onishi, T, Abe, R., et al. 1999, PASJ, 51, 751 

Heiles, C. 1979, ApJ, 229, 533 

Jones, E., Oliphant, T, Peterson, P., et al. 2001-, SciPy: Open source scientific 
tools for Python 

Kalberia, P M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775 



8 



Richard Wiinsch et al.: The Carina Flare: What can fragments in the wall tell us? 



10" 



10- 



10" 



10 




'Tlpeak= 1000 IVlg 

■4 r i~ T — ^ I I d 



20-17 10"^^ 10"^^ 10"^4 10"^^ 10"^2 10"^^ 10"^° 10"^ 

Pext [dyne cnn-2] 



Fig. 8. Maximum of the mass spectrum, mpeak, (solid isolines) and the typical fragment separation, c/peak, (dashed isolines) given by 
PAGI as a function of the shell surface density E and the external pressure P^^.^ . The thick lines show the maximum of the CMF, 
wiCMF, and typical separation of clumps, cImst, obtained from our '^CO APEX observations. The thin dotted lines show the mcMF - 
t/iviST intersection at E = 9 x 10""* gcm"^ and P^^^ = 1.4 x 10"'^ dynecm"^. 



Koyama, H. & Inutsuka, S.-i. 2002, ApJ, 564, L97 

Ki-amer, C, Stutzki, J., Rohrig, R., & Comeliussen, U. 1998, A&A, 329, 249 
MacLaren, I., Richardson, K. M., & Wolfendale, A. W. 1988, ApJ, 333, 821 
Mai'z Apellaniz, J. & Ubeda, L. 2005, ApJ, 629, 873 
Oliphant, T. 2006, Guide to NumPy (Trelgol Publishing, USA) 
Rohlfs, K. & Wilson, T. L. 1996, Tools of Radio Astronomy, ed. Rohlfs, K. & 
Wilson, T. L. 

Rosolowsky, E. W., Pineda, J. E., Kauflfmann, J., & Goodman, A. A. 2008, ApJ, 
679, 1338 

Ryu, D. & Vishniac, E. T. 1991, ApJ, 368, 411 

Sonnentrucker, R, Welty, D. E., Thorbura, J. A., & York, D. G. 2007, ApJS, 168, 
58 

Vazquez-Semadeni, E., Ryu, D., Passot, T., Gonzalez, R. P., & Gazol, A. 2006, 

ApJ, 643, 245 
Vishniac, E. T. 1983, ApJ, 274, 152 

Visser, R., van Dishoeck, E. E, & Black, J. H. 2009, A&A, 503, 323 

Walter, E, Brinks, E., de Blok, W. J. G., et al. 2008, AJ, 136, 2563 

Whitworth, A. P, Bhattal, A. S., Chapman, S. J., Disney, M. J., & Turner, J. A. 

1994, MNRAS, 268, 291 
WiUiams, J. P, de Geus, E. J., & Blitz, L. 1994, ApJ, 428, 693 
Wilson, T. L. 1999, Reports on Progress in Physics, 62, 143 
Wunsch, R., Dale, J. E., Palous, J., & Whitworth, A. R 2010, MNRAS, 407, 

1963 



Appendix A: The DENDROFIND algorithm: 
searching the clump tree 

The DENDROFIND package can be found at 
http : //galaxy . asu . cas . cz/~ richard/dendrofind 
It consists of several py thon scripts an d uses the librarie s 
NumPy (Ol iphantI l2006h and SciPy (Ijones et all I200l4 . 
Figure I A. II shows a flow-chart of the algorithm, the input 
parameters are shown by Table lA.ll and the essential data 
structures are listed in TablelA.2l 



A.1. Description of the aigohthm 

The four input parameters Nlevels, Npxmin, dTleaf and 
Tcutoff are described in section |3] Setting Nlevels is 
straightforward: it should have a value such that increases do 
not lead to different results. The value of Npxmin depends on 
the resolution of the instrument (or the numerical method) and it 
should correspond to the number of pixels of the smallest resolv- 
able structure. Finding values of dTleaf and Tcutoff is not so 
obvious. Fortunately, the dependence of results on these two pa- 
rameters is quite intuitive. Decreasing the value of dTleaf leads 
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START DENDROFIND 



data input 



Tgrid (3d array) 



clump_grid = 


-1 (3d array) 


Tl = 


Tmax 




TL Tgrid 

clumps, clump__grid 



— Tl, Tgrid, clump_grid 
grow_and_connect_clumps clumps 
— new_clump_grid, changes 



n := Tl - (Tmax - 
TcutolT )/( Nlevels ) 



_t_ 



clump_grid := new_clump_grid 



delete_small_c lumps 



compactify_clumps 



- Npxmin 

■ clumps, clump_grid 



clumps, clump_grid 



data output 



draw dendrogram 



' files: clumps. clump_grid 



functions linkage and dendrogram (SciPy) 




END DENDROFIND 



START grow_and_connect_clumps 




changes := 


changes + 1 






px_neighbour_clumps 




-px, clumps, clump_grid 
- neigh_clumps 



- px, neigh_clumps 

- clumps, new_clump_grid 



meree_clumps 
upda[e_touching 



dTleaf, neigh_clumps 
clumps 

new_cliimp_grid 



""Tl, neigh_cliimps 
clumps 




Fig. A.l. The DENDROFIND flow ctiart. The input parameters are printed in red color. Left and right arrows with variables next to 
the boxes with algorithm functions show input and output to/from this function, respectively. Double head arrows denote variables 
used for both input and output. 



to higher numbers of identified structures (as it is easier for lo- 
cal maxima to become independent clumps); and decreasing the 
value of Tcutoff leads to higher fractions of the datacube as- 
signed to clumps (and in most cases also to higher numbers of 
clumps). In many cases, reasonable values for both parameters 
can be set to some multiple of the noise dispersion, e. g., dTleaf 
= Tcutoff - ScTnoise as donc in this work. 

The most important data structures for the algorithm are the 
two 3D arrays Tgrid and clump_grid and the list clumps. The 
Tgrid array is a copy of the PPV data cube including the bright- 
ness temperature T\, for each pixel px. The clump.grid array 
has the same dimensions as Tgrid, and it includes, for each 
pixel px, the ID number of a clump to which pixel px is as- 
signed (-1 means that the pixel is not assigned to any clump). 
The list clumps consists of all the identified clumps represented 
as data objects that include all necessary information. In partic- 
ular, each object contains the clump ID number, the position and 



value of the clump peak temperature, the list pixels of pixels 
that are assigned to the clump, the list touching of clumps that 
directly touch the clump, and the list connected of clumps that 
are connected (even through other clumps) to the clump. The 
lists touching and connected also record information about 
the temperature at which the clump is touching or connected to 
other clumps. Another important piece of information given is 
the list neigh that includes ID numbers of clumps neighbouring 
the processed pixel. In a recent version of the code, neighbours 
of pixel px are defined as pixels which share with px a face, i.e. 
each pixel has six neighbours. A simple single modification of 
the code, however, can change the behaviour so that any other 
pixels (for instance diagonal ones) are also taken into account as 
neighbours. 

The algorithm starts by reading the data from the fits file 
and creating the array Tgrid (see Figure IATI ). Then the array 
clump_grid is created and the value -1 is assigned to all its 
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input param. 


description 


data structure 


description 


Nlevels 
Npxmin 
dTleaf 
Tcutoff 

Table A.l. Lis 


number of temperature levels 

min. number of pixels of a clump 

min. length (in temperature) of a dendrogram leaf 

min. Tb considered for assignment to clumps 

t of input parameters. 


Tgrid 

clump.grid 

clumps 

neigh 

len(neigh) 

Tl 

px 


3D array (PPV) of the brightness temperature 
3D array (PPV) of clump ID numbers 
list of clump objects 

IDs of clumps neighbouring with a processed px. 
length of the neigh array 
actual brightness temp, level in the main loop 
index of the processed pixel (in the code i,j,k) 



pixels. Next, the maximum brightness temperature Tmax is de- 
termined and the main loop of the algorithm is entered. It exe- 
cutes Nlevels iterations decreasing the variable Tl from Tmax 
to Tcutoff. 

In each iteration of the main loop, all pixels px with 
Tgrid(px) > Tl and clump_grid(px) = - 1 are processed in 
several sweeps through the Tgrid array. In the first sweep, a 
new clump is created for each processed pixel which is a lo- 
cal maximum (i. e., Tgrid(px) is greater than Tgrid of all 
px neighbours). An unused ID number is found and assigned to 
clump_grid(px) and a new clump object is created and added 
to the list clumps. In subsequent sweeps, processed pixels which 
are not local maxima are assigned to existing clumps and con- 
nections among clumps are found (see description of function 
grow_and_connect_cltamps below). To make the procedure in- 
dependent of the order in which the pixels are processed in one 
sweep, information about pixel assignments to existing clumps is 
written into a copy of the clurap_grid array new_clump_grid. 
At the end of the sweep, new_clump_grid is copied back to 
cltamp_grid and another sweep is made until no assignment oc- 
curs (this is controlled by the variable changes). 

The grow_and_connect_clumps function consists of sev- 
eral steps executed for each processed pixel (i. e., one with 
Tgrid (px) > Tl and clump_grid(px) = -1). At first, all 
clumps to which neighbour pixels of the processed pixel be- 
long are identified and stored in the list neigh. If neigh in- 
cludes at least one clump, the processed pixel px is assigned to 
the clump from neigh whose peak is closest to px. Further, if 
the neigh list includes more than one clump, four steps are ex- 
ecuted. Firstly, clumps from the neigh list with Tpeak - Tl 
< dTleaf, where Tpeak is the peak temperature of the clump, 
are merged with the clump from the neigh list with the high- 
est Tpeak. Then, the lists touching of all clumps in the neigh 
list are updated so that each clump includes all other clumps 
from the neigh list on its respective touching list. Next, the 
interconnected list is found of all clumps that can be reached 
from the pixel px by moving only through neighbour pixels with 
temperatures greater than Tl. Finally, the connected lists of all 
clumps are updated so that each clump includes all the other 
clumps from the interconnected list on its connected list. 

At the end, after the Tcutoff limit is reached, clumps 
with fewer than Npxmin pixels are deleted and all clumps are 
renumbered so that their IDs are consecutive. Then, the array 
clump_grid and the list clumps are written into files. Finally, 
the SciPy functions linkage and dendrogram are used to cre- 
ate the dendrogram. The function linkage organizes clumps 
into hierarchical clusters based on their connected lists. The 
function dendrogram actually plots the dendrogram finding the 
exact positions of the dendrogram lines in the figure. 

A.2. Sensitivity to noise 

To estimate how sensitive the DENDROFIND algorithm is to 
the level of noise present in the data, we calculate the '^CO 
clump mass spectrum for three different levels of noise (see 



Table A.2. List of variables important for the description of the 
algorithm. 




Fig. A.l. Dependence of the clump mass spectrum on the noise 
in the datacube of the brightness temperature. The solid line 
shows the original mass spectrum (with ctj- = 0.19 K), the 
dashed and dotted lines show spectra obtained from datacubes 
where the noise was artificially increased by factor 2 and 3, re- 
spectively. The dash-dotted line shows the PAGI mass spectrum 
for the same parameters as in Figure|5] 



Figure IA.2I ). The solid line shows the original CMF with the 
noise level ctt - 0. 19 K, the dashed and dotted lines show CMFs 
obtained for ctj- = 0.38 K and ctt - 0.57 K. The resulting spec- 
tra differ mainly in the number of objects (higher ctj leads to 
lower numbers of identified clumps), but the slope of the spec- 
trum and its maximum (mpeak) remain approximately the same. 
There are two reasons for the dependence of the CMF on the 
noise level. Firstly, the noise modifies the brightness tempera- 
ture datacube and results in slightly different masses, positions 
and interconnecting brightness temperatures of clumps. The sec- 
ond, more important reason is that the noise level determines 
values of parameters Tcutoff and dTleaf which are both set 
to 30-7-. We conclude that the CMF slope and peak mass are not 
significantly changed if the noise is increased by up to a factor 
of 3. 



A.3. Sensitivity to ttie velocity resolution 

We also test how sensitive the DENDROFIND algorithm is to 
the velocity resolution. We do it by comparing the CMF ob- 
tained from the original data to CMFs resulting from datacubes 
where two or three velocity channels were added together In 
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Fig. A.3. Dependence of the clump mass spectrum on the veloc- 
ity resolution in the datacube. The solid line shows the original 
mass spectrum (dv - 0.66 kms"'), the dashed and dotted lines 
show spectra obtained if 2 and 3 velocity channels, respectively, 
are added together The dash-dotted line shows the PAGI mass 
spectrum for the same parameters as in Figure |5] 

Figure |A3J the solid line shows the original CMF with the ve- 
locity channel width dv - 0.66 kms ', the dashed and dotted 
lines show CMFs from datasets in which the velocity resolution 
was reduced by factor of 2 and 3, respectively. The CMF with 
twice reduced velocity resolution still shows high mass end that 
is consistent with a power law of slope 2, and a peak at ap- 
proximately ~ 20 Mq. The CMF with dv - 2 kms"', however, 
does not agree with a power-law distribution at all. 

The reason for this behaviour is that the velocity resolution 
of our APEX observation is already quite low, and most clumps 
have FWHM,. smaller than dv - 1.34 kms"' (the original ve- 
locity resolution reduced by 2). Therefore, degrading the veloc- 
ity resolution further results in a drop of the noise in the dat- 
acube (and consequently decrease of parameters Tcutoff and 
dTleaf). The signal-to-noise ratio, however, is not improved. 
Therefore, as the velocity resolution decreases, the algorithm 
tends to find more and more artificial structures. 

Since most of clumps span only a few velocity channels, and 
it rarely happens that the clumps projected onto the /^j-plane 
overlap (see Figure |6] right panel), it is quite improbable that 
one real object is misidentified as several objects due to inter- 
nal velocity flows. We cannot exclude the possibility, however, 
that some clumps could split into several less massive objects at 
higher velocity resolution. 
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Appendix B: Clump properties and individual 
velocity channel maps 

In this appendix we present full versions of Tables [1] and |3] and 
individual velocity channel maps together with maps showing 
identified clumps (see Figures IBTMB.SI i. 
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Table B.l. CO(J=2-l) clump properties - lower region. The meaning of the columns and the means by which the values were 
calculated are explained in |4] Clumps marked by the asterisk at their label are excluded from the CMF, because they lie at the 
border of the observed area. 
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285:59:25.1 


5:20:13.7 


-23.41 


12.9 


0.24 


92.2 


129.9 


0.140 


16E6* 


0.32 


0.48 


0.75 


0.3 


2.3 


285:56:44.8 


5:11:26.8 


-23.41 


3.0 


0.13 


34.5 


60.9 


0.087 


16E7 


0.52 


1.10 


1.02 


0.7 


2.0 


286:02:28.5 


5:11:39.4 


-24.08 


15.8 


0.34 


132.9 


173.2 


0.119 


16E8 


0.34 


0.97 


1.00 


0.5 


1.8 


286:02:57.9 


5:13:47.5 


-23.41 


7.0 


0.23 


90.9 


120.4 


0.076 


16E9 


0.64 


0.34 


0.75 


0.4 


1.5 


286:03:41.9 


5:16:38.3 


-23.41 


3.0 


0.15 


40.1 


71.5 


0.074 


16F1 


0.61 


1.06 


1.03 


0.8 


3.2 


286:01:01.6 


5:08:48.7 


-23.41 


23.8 


0.38 


153.5 


197.6 


0.155 


16F2 


0.81 


1.63 


1.04 


0.9 


2.8 


286:02:13.3 


5:09:02.6 


-24.08 


31.1 


0.44 


181.6 


232.0 


0.171 


16F3* 


0.46 


0.42 


0.85 


0.4 


1.9 


286:00:17.5 


5:05:00.6 


-24.08 


2.8 


0.15 


48.3 


75.4 


0.057 



Table B.2. CO(J=2-l) clump properties - upper region. The meaning of the columns and the means by which the values were 
calculated are explained in ^ Clumps marked by the asterisk at their label are excluded from the CMF, because they lie at the 
border of the observed area. 




Fig. B.l. Clumps identified in the lower region (left two panels) and corresponding brightness temperature (right two panels) in the 
first three velocity channels. Filled regions show clumps 1-18, dotted regions show clumps 19-35. Colors of clumps correspond to 
colors of dendrogram leaves in Figure|2l 



Richard Wiinsch et al.: The Carina Flare: What can fragments in the wall tell us?, Online Material p 4 




Fig. B.2. Clumps identified in the lower region (left two panels) and corresponding brightness temperature (right two panels) in the 
second three velocity channels. Filled regions show clumps 1-18, dotted regions show clumps 19-35. Colors of clumps correspond 
to colors of dendrogram leaves in Figure |2] 




Fig. B.3. Clumps identified in the lower region (left two panels) and corresponding brightness temperature (right two panels) in the 
last three velocity channels. Filled regions show clumps 1-18, dotted regions show clumps 19-35. Colors of clumps correspond to 
colors of dendrogram leaves in Figure|2l 




Fig. B.4. Clumps identified in the upper region (left two panels) and corresponding brightness temperature (right two panels) in the 
first three velocity channels. Colors of clumps correspond to colors of dendrogram leaves in Figure [3] 



Richard Wiinsch et al.: The Carina Flare: What can fragments in the wall tell us?, Online Material p 7 




Fig. B.5. Clumps identified in the upper region (left two panels) and corresponding brightness temperature (right two panels) in the 
last two velocity channels. Colors of clumps correspond to colors of dendrogram leaves in Figure|3] 



